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Abstract 

The nonleptonic two body decays D — > PP and D — ► PV are investigated based on the 
diagrammatic decomposition in a generalized factorization formalism. It is shown that to fit the 
experimental data, the SU(3) flavor symmetry breaking effects of the coefficients a^s should be 
considered in D —* PP decay modes. In D — > PV decays, the final state hadron structure due 
to the pseudoscalar and vector mesons has more important effects on the coefficients a^s than 
the SU(3) symmetry breaking effects. It is found that the nonfactorizable contributions as well 
as that of the exchange and annihilation diagrams are significant in these decays. 
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I. INTRODUCTION 



Charmed meson nonleptonic two body decays have been an interesting subject of re- 
search 0, 0, 1110] f° r a long time as it can provide useful information on flavor mixing, 
CP violation |5[ and strong interactions. The theoretical settlement of this transition 
type generally appeals to the factorization hypothesis. Empirically, nonfactorizable cor- 
rections which result from spectator interactions, final state interactions and resonance 
effects should be considered. The nonfactorizable corrections are believed to be significant 
, and they are relatively hard to be calculated because the charmed quark is not heavy 
enough to apply the QCD factorization approach [tJ or PQCD approach in a reliable 
manner. Fortunately, a great number of precise experimental data on charmed meson 
nonleptonic two body decays have been accumulated in recent years. Many new results 
are expected soon from the dedicated experiments conducted at BES, CLEO, E791, FO- 
CUS, SELEX and the two B factories BaBar and Belle. Phenomenological models based 
on all kinds of symmetries are of quite importance to guide the theoretical studies and 
explore new physics (HEll II |- But in some cases, the symmetry breaking effects can be 



significantly enhanced. 

In the quark diagrammatic scenario, all two-body nonleptonic weak decays of charmed 
mesons can be expressed in terms of six distinct quark-graph contributions (H[I2 : U) a 
color-favored tree amplitude T, (2) a color-suppressed tree amplitude C, (3) a W-exchange 
amplitude E, (4) a W-annihilation amplitude A, (5) a horizontal W-loop amplitude P 
and (6) a vertical W-loop amplitude D. The P and D diagrams play little role in practice 
because the CKM matrix elements have the relation V* s V us ~ —V* d V u d which will result 
in cancellations among these diagrams. 

Based on SU(3) flavor symmetry, the T, C, E and A amplitudes were fitted from the 
measured D meson decay modes jiy, EH • These amplitudes help one to understand the 
generality of charmed meson decays. But since SU(3) flavor symmetry breaking effects 
appear to be important |l3|, E3, these fitted data can not describe the specific properties 
in certain decay modes. In 15], we investigated in detail both the Cabibbo-allowed and 
singly Cabibbo-suppressed D — > PV decays based on the diagrammatic decomposition 
in the factorization formalism and found that the SU(3) symmetry breaking effects in 
the D — > PV decays are significant. Two sets of solutions were found in the formalism 
of factorization. The case (I) solution can provide satisfactory explanation in a natural 

manner on the process D + — ► K°K* + which is thought to be a puzzle ^3- -^ u ^ the 
solution is hard to be explained from the theoretical point of view because this solution 
requires such an unexpected large correction from nonfactorizable contributions that the 
strong phase of a^ p has around 150° deviation from that of Wilson coefficients c%. The 
case (II) solution shows relatively small correction from nonfactorizable contributions and 
hence seems more reasonable from theoretical point of view. But the solution predicts 
a relatively small branching ratio of the process D + -> K K* + in comparison with the 
experimental result. With such a treatment via solving fifteen equations for extracting 
out the same numbers of parameters, it is hard to consider the experimental uncertainties 
in To investigate what impacts the experimental uncertainties will bring to the 

extracted parameters, it is useful to make a systematic analysis with taking into account 
the experimental uncertainties. 

In this paper, we will perform a \ 2 fitting procedure on charmed mesons decaying to 
a pseudoscalar and a vector meson (D — > PV) and also decaying into two light pseu- 
doscalar mesons (D — > PP) by using the quark-graph description based on a generalized 
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factorization formalism which reflects SU(3) flavor symmetry breaking effects. Firstly by 
dividing these diagrams into factors including SU(3) flavor symmetry breaking effects and 
introducing parameters describing the overall properties, we arrive at two sets of solutions 
for the parameters from fitting experimental data. In the view point of diagrammatic de- 
composition, the generalized QCD parameters ai(i = 1,2) will be classified into two sets 
of parameters af and aj . The difference between af and a( arises from the final state 
hadron structure of the pseudoscalar and vector mesons in D — > PV. In D — > PP decays, 
we will show that, to fit the experimental data, one should classify the parameters into 
af and af, which means that the SU(3) flavor symmetry breaking effects are important and 
need to be considered in the a^s. Thus we can arrive at a conclusion that the coefficients 
a\ and a-i depend on either the final state hadron structure or SU(3) flavor symmetry 
breaking effects. For D — > PP decay modes, the SU(3) flavor symmetry breaking effects 
play an important role in the coefficients a\ and a 2 , while for D — > PV decays, the final 
state hadron structure becomes more important for the contributions to the coefficients 
than the SU(3) symmetry breaking effect does. The contributions of SU(3) flavor sym- 
metry breaking effects to ai and a 2 can be neglected in D — > PV decay modes. Using 
the fitted parameters as inputs, we are led to predictions for the branching ratios of other 
decay modes which are expected to be measured in the future. In studying the breaking 
of the SU(3) symmetry relations, we are able to quantify the SU(3) breaking effects. The 
breaking amount of the SU(3) symmetry relations in some channels can be significant so 
that it becomes unreliable to use the SU(3) relations to make predictions for some decay 
modes. 

The paper is organized as follows. In section m we list the flavor decomposition of the 
corresponding mesons and present the quark-diagram description for the relevant decay 
modes. In section llll[ the parameterized formalism based on factorization is introduced 
to investigate the processes. We then perform a fit procedure in section II VI to extract 
the parameters and present predictions for thirty three D — > PP decay modes and sixty 
two D — > PV decay modes. The SU(3) flavor symmetry breaking effects are discussed in 
section A short summary and remark is given in the last section. 



II. NOTATION AND QUARK-DIAGRAM FORMALISM 



We adopt the following quark contents and phase conventions which have been widely 
used HI H Hill. 

• Charmed mesons: D° = —cu, D + = cd, = cs; 

• Pseudoscalar mesons P: tt + = ud, tt° = (dd — uu)/y/2, tt~ = —du, 
K + = us, K° = ds, K° = sd, K = —sit, 77 = (—uu — dd + ss)/y/3, 
rj = (uu + dd+ 2ss) / v6; 

• Vector mesons V: p + = ud, p° = (dd — uu) / y/2, p~ = —du, u = 
(uu-\- dd)/y/2, K* + = us, K*° = ds, K* {) = sd, K* = — su, (f) = ss. 

In the above notations, u, d and s quarks transform as a triplet of flavor SU(3) group, 
and — u, d and s as an antitriplet, so that mesons form isospin multiplets without extra 
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signs. In general, the rjr]' mixing are defined as 





with 7] = (uu + dd + ss)/y3 and r]s = (—uu — dd + 2ss)/v / 6- For convenience, we 
have taken the mixing parameter as 6 = 19.5° = sin _1 (l/3) which is close to the value 
= 15.4° extracted from experiment 48j. 

The partial width F for D — > PP and P — > PV decays is expressed in terms of an 
invariant amplitude A. One has 



for D -> PP and 



for P — > PV, where 



P 



,3 



r(p^PV) = -^- T |^| 2 (3) 



8ttM 



P 



(M|, - (mj + m 2 ) 2 )(M 2 D - (m x - m 2 ) 2 ) 



2M, 



denotes the center-of-mass 3-momentum of each final particle. 

In P — > PP decays, to describe the flavor SU(3) breaking effects, a subscript s or d is 
attributed on T and C diagrams to distinguish the initial c quark transits to s quark or d 
quark. The subscript s or d is attached to the diagrams E and A dominated by the weak 
process cq 1 — ► g2<?3 when the antiquark g 3 is s or d. In P — » PK decays, a subscript P 
or V is assigned to T and C, which are induced by c — > qzQ\<i2 with the spectator quark 
containing in pseudoscalar or vector final meson. The subscript P or V is labelled to E and 
A graphs which are dominated by the weak process cq 1 —>■ q2q 3 when the final antiquark 
q 3 stays in the pseudoscalar or vector meson. S is added before E or A to distinguish the 
exchange or annihilation graph involving in final singlet state contributions which result 
from disconnected graphs. The total contributions of the SE and SA graphs involving 
in 7t° and p° mesons are equal to zero because their contributions resulting from uu and 
—dd offset each other due to the isospin SU(2) symmetry. In the numerical analysis, 
we will assume that the contributions of the SEp and SEy graphs involving in u and 
mesons are negligibly small since they seem not to contradict with the Okubo-Zweig- 
Iizuka rule. But the amplitude SA V seems to play an important role in the P+ — > p + i] 
and P+ — > p + r]' processes [l9[. In the ideal mixing case, the process P+ — > n + u! has the 
amplitude representation as -^(Ay + Ap + 2SAp). Since u> has a similar quark structure 
in comparison with rj and rj', we assume that SAp has an important contribution in 
P+ — > 7i + u. In the present paper, we shall not consider the processes which receive 
contributions from SAy and SAp diagrams resulting from the final state particles rj, 7/ 
or uu. The sign flips in the presentation of some relevant Cabibbo-favored modes, as well 
as that of some doubly Cabibbo-suppressed modes, come from the quark contents of final 
light mesons. In the singly Cabibbo-suppressed modes, the sign flips may come either 
from the quark contents of the final light mesons or from the CKM matrix element V* d V u d 
since V* S V US ~ —V* d V u d and we choose V* S V US in the calculations. In Table HVI and Table 
rVT a prime and double prime are added to the diagrams of singly Cabibbo-suppressed 
modes and doubly Cabibbo-suppressed modes respectively to distinguish them from the 
Cabibbo-favored ones. 
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III. FLAVOR SU(3) SYMMETRY BREAKING DESCRIPTION IN GENERAL- 
IZED FACTORIZATION FORMALISM 

To investigate the SU(3) flavor symmetry breaking effects, we take the formalism of a 
generalized factorization approach 2, 20]. 

For D — > PP decays, amplitudes can be written in the form as 

Ts, d = ^V qiq2 V; qs a Ta J Pl (ml - m\)F^{rn\) : (4) 

C s , d = ^V qiq2 V^ac a J Pl (m 2 Di - m 2 p 2 )F D ^(m 2 Pi ), (5) 
G F 

E s,d = -^ V qiq3 V cq 2 a Es,dfDii ( 6 ) 
Gp 

A s,d = —^VqiqaV^A^jDi- (7) 

For D — > PV decays, amplitudes can be written in the form as 

Tv = ^V qiq2 V; q3 a Tv 2f P m Di A^ v (m 2 P ), (8) 

Tp = ^V qiq2 V: q3 a Tp 2f v m Dt F^ p {m 2 v ), (9) 

Cv = |^K;a Cf 2K^K), (10) 

C P = ^V qiq2 V: qs a Cp 2f v m Di F^ p {m 2 v ), (11) 



E 



G 

= -|^i?3K>^,p 2 /A m A» ( 12 ) 



V,P 



V2 



Gp 

a v,p = -j^V q2q3 V* qi a Avp 2f D .m D .. (13) 

Di denotes D ± , Dq or D s . Fq, F\ and Aq are formfactors defined in the following formalism 

(P(p)\ffc\D(p D )) = \( PD +p), - m2D ~ m2p A F.iq 2 ) + ^'^ q^iq 2 ), (14) 

(V(p)\qY(l-l 5 )c\D(p D )) = - t (m D + m v )A 1 (q 2 )(e^ - ^V) 

qz 

+j ^!L (£ . . qK(PD + rY _ _ . )ilohV 

mo + mv Q Q 

- 2V{ / ] e^e* aPDp p 1: (15) 
m D + m v 

with q = pp, — p. fp and fy are decay constants defined as 

(P(p)|gi7 M 75g 2 |0) = -ifpjf, (16) 
(V(p)\q 1 ^q 2 \Q) = f v m v e». (17) 



In naive factorization hypothesis, one has the following equalities 



ot s = a Td = a Tv = a Tp = ai(fx), (18) 
«c s = «c d = a Cv = a Cp = a 2 (/x), (19) 



with 



oi(^) = C M + ^ca(^)» (2°) 

02^) = c 2 (/i) + ^Ci(/i), (21) 

denoting the relations between quantities ai j2 and Wilson coefficients ci )2 . iV c is the 
number of colors, /i is the renormalization scale at which c\ and c 2 are evaluated. So ai 
and a 2 are common real quantities of a certain process in quark level. To be more explicit, 
for decay modes induced by c — > s transition, a\ and a 2 are invariant among all modes in 
naive factorization hypothesis. 

However, naive factorization approach meets difficulties in describing all charmed me- 
son decays, particularly for the decay modes which involve in the color-suppressed dia- 
grams due to the smallness of |a 2 |. Furthermore, the coefficients a\ and a 2 in Eqs. ()18|) 
and (|19|) depend on the renormalization scale and 75 scheme at the next to leading order 
expansion. It is necessary to consider the nonfactorizable corrections which involve in 
hard spectator interactions, final state interactions and resonance effects etc. In a general 
case, one can express a\ and a 2 in the form 

ai(p) = ci(p) + (— +XiO))c 2 (/V), (22) 

a 2 (» = c 2 (/x) + (-^r + X2(/i))ci(/i), (23) 

with Xi(/-0 an d X^in) terms denoting the nonfactorizable effects. With nonfactorization 
corrections the equalities (fTSj) and (fTU|) are not yet satisfied because each a, should contain 
terms from different corrections. The nonfactorization corrections can also bring phase 
differences among these coefficients, and then ajS (i = T s>d> y ;P , C St d,v,p) become complexes. 
Currently, explicit calculations of total nonfactorizable corrections are not yet possible. 
In D — > PV decays, we shall take all OjS as independent complex parameters and assume 
that the corrections do not depend on individual decay process at certain scale. In other 
words, we do not consider SU(3) flavor symmetry violation contributions to <2jS and it 
is supposed that mass factors, decay constants and formfactors have taken on the whole 
SU(3) symmetry breaking effects. While in D — > PP decays, the mass factors, the form 
factors and decay constants fail to account for the large SU(3) flavor symmetry breaking 
effects in D — > tttt, D — > ttK and D — ► KK. Nonfactorizable contributions may cause 
large SU(3) symmetry breaking Two sets of coefficients and af are introduced to 



describe the SU(3) flavor symmetry breaking effects induced by nonfactorizable contribu- 
tions. In both D — > PP and D — > PV decays, the SU(3) symmetry breaking effects are 
not considered in the strong phases in our present analysis. 

The exchange and annihilation diagrams have the following expressions in naive fac- 
torization approach: 

E s , d = ^V qiq3 V^a n J s J Di (m 2 Pl - m 2 P2 )F p ^(m 2 Di ), (24) 
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A s , d = ^V^V^at joXm^ - m%) F™ (m 2 Di ) , (25) 
G f 



E 



v, 



p = ^V qiq3 V; q2 a% p 2f Di m D A^(m 2 D J, (26) 
A v ,p = ^V q2g3 V^y/ Vp 2f DimDi A^(m 2 Di ). (27) 

The formfactors FQ lBl {m 2 D ) and Aq V (m? Di ) involving in the above formula are not mani- 
fest to relate directly to experimental measurements. The factorizable contributions of the 
exchange and annihilation diagrams are believed to be small. The main contributions of 
these diagrams may result from the nonfactorizable forms. Through intermediate states, 
these diagrams relate to the tree diagram T and color-suppressed diagram C jU 2^ 



Their contributions may be important and can not be ignored. In our present considera- 
tions, we use dE^Aiii — s > d, V, P) defined in eqs.flOJ), ©, CH) an d (fT^j) as global parameters 
to describe mainly the nonfactorizable contributions. By these definitions, the parame- 
ters ciEi,Ai w iU have two dimensions of energy in D — ► PP and will be dimensionless in 
D -> PV. 



IV. NUMERICAL ANALYSIS AND RESULTS 



The explicit evaluation of the relevant formfactors in the factorization formula P|). 
© and (J%j)- (fTT]) is a hard task because of the nonperturbative long distance effects of 
QCD. Various methods, such asQCD sum rulesji^, 24]|, lattice simulations |25|, |26| and 
phenomenological quark model 27], 28[, have been developed to estimate the long distance 
effects to rather high certainties. The formfactors of D mesons decaying to light mesons 



have been widely discussed in [29, 30, |3JJ, |32|, |33|, |34]. In our present considerations, we 
shall use the results of form factors obtained by Bauer, Stech and Wirbel 0, based on 
the quark model. They have been found to be rather successful in describing a number 
of processes concerning D mesons. The values of the relevant formfactors evaluated at 
q 2 = are listed in Table |l] For the dependence on q 2 , the formfactors are assumed to 
behave as a monopole dominance 



D - P : F (q 2 ) = (28) 

1 — q z /m F «* 

Fl(q 2 ) = : Fl if 2 , (29) 
1 — q l jm F , 

D^V: A (q 2 ) = - Ao ^ 2 , (30) 



where rap, m F * and m F ** are the pole masses given in Table 

It is noted that the formfactors are more appropriate to be viewed as the relative scaling 
factors that characterize one source of SU(3) flavor symmetry breaking effects in hadronic 
matrix elements since we take the a^s as free parameters that need to be extracted from 
experimental inputs in the present method. The relative ratio between the formfactors is 
what we really care about. 

The input values for the light pseudoscalar and vector decay constants are presented in 
Table HP pil l36l| . These values generally coincide with experiments. The decay constants 
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/* f™, and /*/ involving in factorization formula should be defined as follow |35[ 

(o\uYi 5 u\^'\p)) = zf; (f) p^ (31) 

(0\srj 5 s\rt , \p))=if^ )P >*. (32) 
Then the quantities fL /"/ and /*/ take the formalism 

^ = vl COS( ^ + vl sin< ^' ^ 

/„ = ^= cos - sin 0, (34) 

Making use of these definitions, the following factorization formalisms are adopted in the 
D — > r](i]')V transition calculation 

2C V (A -> 17V) = ^§^ 1?2 ^ 3 a Cv 2(/« + f°)m Di A^ v (ml), (37) 
72 



CV(A -> t/V) = ^^ 1?2 V c ; 3 a Cv 2(^ - /> D ir F K). (3* 



The other parameters used in the numerical calculation are the masses of relevant 
mesons, lifetimes of charmed mesons and relevant CKM matrix elements. We adopt the 
relevant values given in |37| . 

For convenience, we may express the complex parameters a« as 

a i =\a i \e i5 °i. (39) 

The S ai s characterize the strong phases. One can always choose 5 aTg = in D — > PP and 
$a Tv = in D —> PV so that all the other strong phases are relative to 5 aTs and 5 aTv . 
There are 15 independent parameters to be extracted from experiments in both D — > PP 
and D -> PV. 

To conduct a fit procedure, we construct a x 2 function which has the following form 

x» = £ ifAai) ~, U ' )? (40) 

3 a 3 

where (fj) and <jj are the central values and corresponding errors of the experimentally 
measured observables. fj(di) are the theoretical expressions for the observables. They are 
the functions of parameters djS. The set of a^s which minimizes the x 2 function will be 
regarded as the best estimated values. 

There are 17 experimental data points for 13 parameters (a s A does not appear in these 
experimental data) and 22 data points for 15 parameters, as shown in Table HVI and Table 
13 respectively. We list in Table ITTT1 the parameters with la errors obtained in our present 
analysis. FIT a and FIT A are obtained without any constraint to the parameters. A 
large |a|/af| ~ 2.0 ratio predicted by FIT a is an indiction of inscrutably large flavor 
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SU(3) breaking effects. Constraining the ratio to the smallest extent, we get FIT (3 
with the ratio lal/ ^! ~ l-l- By "the smallest extent", we mean that, if we continue to 
suppress the ratio down, the predicted branching ratios of some decay modes in Table 
IIVI will be inconsistent with the experimental data. FIT A predicts an unusually large 
ratio (02/01 | ~ 1.1 which indicates that the nonfactorizable contributions to a\ are of 
great importance. By constraining the value of \aX\ to be as small as possible, we obtain 
FIT B with the ratio 0.9. The next leading order Wilson coefficients ci(m c ) = 1.174 and 
C2(^ c ) — —0.356 in the naive dimensional regularization (NDR) scheme or Ci (m r ) = 1.216 
and C2(m c ) = —0.424 in the 'tHooft-Veltman (HV) scheme are given in Ref.j38( when 
AjjTg = 0.215GeV. The present relatively large values of | \\ and \a^\ can not 

be explained from formula (|2Tj) . They imply that nonfactorizable contributions are of 
significance in both D — > PP and D — > PV decays. To fit the experimental result of 
Br(D° -»■ K°K°) = (0.071 ± 0.019)%, a% should differ much from a%. In D -»• PV 
decays, because we have considered the errors of experimental data in the \ 2 fit and used 
more experimental results as constraints, the present resulting parameters appear more 
reasonable than that of case (I) solution presented in [Hj], as the strong phases of the 
parameters af' y and are not in contradiction to that predicted from QCD. 

We present the resultant predictions for a variety of charmed meson decay processes 
in Table HVl for D — > PP decays and in Table IVl for D — > PV decays. Note that there are 
no enough experimental data to extract the parameter a s A in D — > PP decays. To make 
predictions for some relevant decay modes which receive contribution from A s diagram, 
we take the assumption a s A = a d A . For a comparison, we also list the results obtained in 
Ref . [ioT] . The predictions for a number of singly and doubly Cabibbo-suppressed modes 
can be used to test our present analysis in the near future. 

Note that in the assumption of SAp = 0, we have the branching ratio 14% for the 
process — > tt + lj, which is much larger than the experimental result (0.28 ± 0.11)%. 
To accommodate the experimental data, significant contribution from SAp diagram, i.e. 
SAp ~ — A p , should be introduced [l5| . 



V. SU(3) FLAVOR SYMMETRY BREAKING 



As pointed out in Ref . |l3j, |l4| , SU(3) breaking effects in charmed meson decays appear 
to be important. The violation may come from the finite strange quark mass, the final 
state interactions and resonances. In the SU(3) flavor symmetry limit, there are a number 
of relations among different decay modes. Based on the above extracted values for the 
parameters, we can discuss how large are the SU(3) breaking effects in D — ► PP and 
D — > PV decays. 

We present these relations in Table IV1I for D — > PP and Table IV11I for D — > PV. 
The left hand side(LHS) values of the relations whose deviation from unit represents the 
breaking amounts of SU(3) flavor symmetry relations are listed in the second columns. 

It is noted that though these relations deviating from unit reflects the SU(3) flavor 
symmetry breaking effects, the ones composed of three decay modes and those composed 
of two decay modes have different sources of breaking terms. To be clear, we take the 

expressions ^^^(li^X)^^ and ^kP^W as examples. 

\XA(D + -> Tr+iT ) + y/2A(D + -> 7r+p°)| 
\\V2A(D+ ^n+p°)\ 
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\{Ty + C P ){D + -> 7T + F*°) - (7V + Cp - + Ay)(£) + -> 7T + p°)| 

|^l(L) -> tf+K"*-)| |Ty(£>° K+K*-) + E P (D° -> 



(41) 
(42) 





-> 71 K. J 




-> 71 K. ) 


A y (L>+ - 


- 7T V) = 


Ap(D + - 


- 7T + P°) = 


Ty(D° - 




E P (D° - 


-> ^IT" 



|,A(D ->7r + /3-)| |-Ty(D°^7r+p-)- J Ep( J D ^7r+p-)| ' 
In the limit of SU(3) flavor symmetry, we have the following relations 

= r v p+- > 7rV), (43) 
= Cp(^ + - ^V), (44) 
A v {D+^n+p°), (45) 
Ap(£>+ - tt+p ), (46) 
= T v (£>°-.7rV), (47) 
= £ P (£> -7rV), (48) 
which make the ratios eq.(|4*T|) and eq.(jHJ) equal to one. But from formula (|Hj)- (fT3^) . 
one can find that relations in eqs. (|43|) - (j4*%|) are in general not valid. The different 
masses of the charmed mesons and the final light mesons, and the different values of 
formfactors and decay constants can break the relations in eqs. (|43j ) -([48| ) . and thus break 
the SU(3) flavor symmetry relations in (|41|) and (|42|) . In addition, by comparing with 
(|41|) and (|42|) . one can see that the relations concerning only two decay modes represent 
the relative SU(3) flavor symmetry breaking amounts of the same diagrams which we call 
the main diagrams for convenience in later use, while the relations consisting of three 
decay modes contain additional SU(3) flavor symmetry breaking effects from the other 
diagrams. So in the relations containing three decay modes, if the SU(3) flavor symmetry 
breaking contributions of the other diagrams have comparable amounts in comparison 
with the main diagrams, then the relations will be broken down badly. The main diagrams 
\T + C\ in D + -> tt + K°, \A V - A P \ in D+ -> n + p° and \T V + C P \ in D + -> 7T + K*° are 
relatively small, which usually leads to a significant breaking for the relations when taking 
A(D + — > tt + T{ ), A(Df — > 7r + p°) and A(D + — > tt + K* ) as denominators. We present 
explicitly some of these relations calculated from the parameters of FIT a and FIT A as 
follows: 

\V2*A(D+ -> K+^) - kA{D+ - K°tt+)\ = 

0.60, (50) 



\XA(D+ - 

\XA(D + -> n+K* ) + y/2A(D + -> 7r+p°)| 
\XA(D+ -> ZV*+) + >/2X(I>+ -> K*+7r°)| 



|AVlA(£>+^7r+p )l 



|A.A(D+ -> jr+K* ) + .A(£)+ -> Ar+Ar* ) 
\\A(D+ ^ 7i+K*°)\ 



1.42, (51) 



0.58. (52) 



It is obvious that SU(3) flavor symmetry analysis is not applicable to such processes. 

Besides the mass factors, the form factors and decay constants, one should also consider 
the contributions of dj factors when studying the SU(3) symmetry breaking effects in 
D — ► PP decay modes. The situations are more complicated than that in D — > PV 
decay modes. General speaking, the SU(3) flavor symmetry breaking effects are more 
important in D — > PP decays. The first two relations in Table IVII and Table I VI II still 
conserve because all the decay modes in them form an isospin triangle respectively. 
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VI. SUMMARY AND CONCLUSION 



We have performed a x 2 fitting analysis on the D — > PP and £) — > decays 
in the formalism of the factorization hypotheses. To fit the experimental data, it is 
vital to consider the SU(3) flavor symmetry breaking effects of the coefficients a^s in 
D — > PP decay modes. In D — > PV decays, the final state hadron structure of the 
pseudoscalar and vector mesons has more important impact on the coefficients a^s than 
the SU(3) symmetry breaking effects. The nonfactorizable contributions, as well as that 
of the exchange and annihilation diagrams, are found to be important in these decays. 
In the formalism of the relations obtained in the SU(3) symmetry limit, the total SU(3) 
symmetry breaking amount of certain processes in D —>■ PP can reach 120% when the 
three symmetry breaking effects due to factors, mass factors and due to form factors 
and decay constants become to be coherently added. The total breaking amount of 
some processes in D — ► PV can add up to 50%. The breaking amount of the SU(3) 
symmetry relations in some channels is so significant that it becomes unreliable to use the 
SU(3) relations to make predictions for some decay modes. More precise measurement on 
the process D + — > K°K* + is important for understanding the SU(3) symmetry breaking 
effects and nonfactorizable contributions. As an independent check, it is useful to measure 
the process — > K°p + . The processes D° — ► 7r + p~, D° — > 7i~p + , D° — > 7r°p°, D + — > 
n+LU, D+ -> n°p + , D+ -> K°p+, D+ -> ir°K*+, D+ -> K+w and D+ -> vr°iC + are 
predicted to be at the experimental sensitivity. It is expected to explore the final hadron 
structure and SU(3) flavor symmetry breaking effects in D — > PP and D — > PV decays 
in BES, CLEO-c, BaBar and Belle. 
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TABLE I: Relevant formfactors at zero momentum transfer for D — ► P and D — > V transitions 
and pole masses in BSW model. 



Decay 


D -> TT 


D - p(uj) 


D -> K 


D ^ K* 


D S ->K 


D s -» K* 




D — > rj/r}' 


D s -» 77/77' 


*l 


0.692 




0.762 




0.643 






0.681/0.655 


0.723/0.704 


A 




0.669 




0.733 




0.634 


0.700 






m F (GeV) 




1.87 




1.97 




1.87 


1.97 






rn F ,(GeV) 


2.01 




2.11 




2.01 






2.01 


2.11 


m F ,*(GeV) 


2.47 




2.60 




2.47 






2.47 


2.60 



TABLE II: Values of decay constants in MeV. 



u 


fK 


h 


/o 


/d 


IDs 


/p 


fK* 


/ w 


/</> 


f D . 


/r>; 


134 


158 


168 


157 


200 


234 


210 


214 


195 


233 


230 


275 
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TABLE III: Parameters a,s fitted from experimental data at la errors. The first entry is for 
amplitude and the second entry for the strong phase. a\' ' ' and a^' ' ' denote aT sdVP and 
a c aA , ViP respectively. 



D -> PP 


D -> PV 


X 2 /d.o.f. 


FIT a 
4.06/4 


FIT j3 
8.16/4 




FIT A 

8.22/7 


FIT B 
10.30/7 


af 


1.08 ± 0.04 


1.10 ±0.03 




1.13 ±0.08 


1.10 ±0.07 


af 


1.04 ±0.09 
(8.73 ±7.96)° 


1.09 ± 0.09 
(11.98 ±7.85)° 


af 


1.29 ± 0.04 
(10.04 ± 16.62)° 


1.29 ± 0.04 
(-1.36 ± 13.52)° 


a| 


-0.73 ± 0.04 
(-26.76 ± 1.60)° 


-0.73 ±0.04 
(-26.25 ± 1.55)° 


aX 


-1.19 ±0.06 
(-11.09 ±20.01)° 


-1.00 ±0.05 
(-10.74 ± 10.31)° 


4 


-0.36 ±0.20 
(-53.40 ± 28.65)° 


-0.65 ± 0.06 
(-35.12 ± 14.50)° 


af 


-0.78 ±0.03 
(-21.75 ± 1.38)° 


-0.77 ±0.02 
(-22.15 ±2.23)° 


a s E (GeV 2 ) 


0.24 ±0.11 
(-49.43 ± 20.63)° 


0.25 ±0.11 
(-58.75 ± 19.48)° 


a\ 


0.07 ±0.03 
(-166.87 ±50.96)° 


0.26 ± 0.04 
(-115.50 ± 12.78)° 


a d E (GeV 2 ) 


1.01 ±0.08 
(-120.94 ±2.92)° 


1.01 ±0.08 
(-122.02 ±2.82)° 


a E 


0.51 ±0.02 
(82.82 ±4.01)° 


0.50 ±0.03 
(78.55 ± 5.73)° 


a s A (GeV 2 ) 






a A 


0.52 ±0.05 
(-76.66 ± 25.78)° 


0.53 ±0.05 
(-78.24 ±26.13)° 


a A {GeV 2 ) 


0.43 ± 0.09 
(90.04 ± 13.75)° 


0.46 ± 0.09 
(89.94 ± 14.27)° 


a A 


0.59 ±0.03 
(-76.29 ±41.25)° 


0.60 ± 0.02 
(-77.93 ±31.52)° 
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TABLE IV: Predicted branching ratios for charmed mesons decaying to two pseudoscalar mesons. 
Single prime and double primes are added to the representations to denote the singly Cabibbo- 
suppressed processes and doubly Cabibbo-suppressed processes. C Sl and C S2 , as well as C dl and 
Cd 2 1 result from the exchange of the final mesons. 



Decay Modes 


Representation 


Experimental 


Present B x 10 2 


LP [401 






B x 10~ 2 


FIT a 


FIT 


B x 10" 2 




K-TT+ 


T s + E d 


3.80 ±0.09 


3.79 


3.80 


3.847 




K°7T° 


71 ~ E d) 


2.28 ±0.22 


2.27 


2.24 


1.310 




K rj 
K rj 


1 r 


0.76 ±0.11 


0.80 


0.81 






- i {C s + 3E d ) 


1.87 ±0.28 


1.85 


1.88 








-(T'd + K) 


0.143 ± 0.007 


0.144 


0.144 


0.151 




vr 7r° 


-^(C' d -E' d ) 

V z 


0.084 ± 0.022 


0.078 


0.097 


0.115 




K+K- 


T' s + E' s 


0.412 ±0.014 


0.413 


0.413 


0.424 




K°K° 


E's-E' d 


0.071 ±0.019 


0.069 


0.062 


0.130 


D° 


K+ir- 


-(T d ' + E») 


0.0148 ± 0.0021 


0.0150 


0.0151 


0.033 




rjTt 


±(C' s + C' di -C> d2 -2E' d -SE') 


— 


0.069 


0.068 






rj'n 


$3 ( 2C 's-C' dl +C > d2 +2E' d +4SE') 


— 


0.088 


0.091 






rjrj 


^(2C' S + 2C' d -2E' S + 2E' d + ASE) 


— 


0.011 


0.016 






rjrj' 


-J=(2C' -C -C -C> -iE>„-1B> -7SE) 
V18 1 °2 a l "2 " 




0.026 


0.030 






K°7T 




— 


0.002 


0.005 


0.008 




K°ri 


-^(C>i-E>l + SE») 


— 


0.001 


0.002 






K°rj' 


MC>> + ZE» + ASE») 

V D 


— 


0.0 


0.0 






K U vr+ 


T s + C s 


2.77 ±0.18 


2.76 


2.76 


2.939 




vr+vr 


-^{T' d + C' d ) 


0.25 ±0.07 


0.25 


0.19 


0.185 




r]TT + 


MT' d + C' a + C' d +2A' d + SA') 


0.30 ± 0.06 


0.34 


0.37 






rj'n + 
K+K° 


-^ T, d -Ws+C' d + 2A> d + iSA') 


0.50 ±0.10 


0.45 


0.42 




D+ 


T' s -A' d 


0.58 ±0.06 


0.62 


0.62 


0.764 




K°TT+ 


-ic'^ + K) 


— 


0.012 


0.026 


0.053 




K+7T 




— 


0.021 


0.023 


0.055 




K+rj 






0.011 


0.012 






K + rj' 


--^{Tt + ZA'i + ASA") 




0.005 


0.006 








C s + A d 


3.6 ±1.1 


3.06 


3.13 


4.623 




n + rj 


j= 3 (T s -2A d -SA) 


1.7 ±0.5 


1.05 


1.09 


1.131 




TT + rj' 


^(T s + A d + 2SA) 


3.9 ±1.0 


4.19 


4.43 




Df 


TT+K° 


-(T' d -A' s ) 


< 0.8 


0.24 


0.26 


0.373 




TT°K+ 


-js(C' d + A' a ) 




0.047 


0.090 


0.146 




TjK+ 


j-<T' a + C' a + C' d -SA') 




0.055 


0.040 


0.300 




7j'K+ 


■±{2T> a + 2C' s -C' d + 3A> 3 +ASA>) 




0.090 


0.102 






K+K° 


-W + c%) 




0.014 


0.010 


0.012 
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TABLE V: Predicted branching ratios for charmed mesons decaying to one pseudoscalar and 
one vector meson. Single prime and double primes are added to the representations to denote 
the singly Cabibbo-suppressed processes and doubly Cabibbo-suppressed processes. 



Decay Modes 


Representation 


Experimental 


Present ^(xlO^ 2 ) 


List I4U| 








I'll A 


T7 1 T r P "D 
Ml D 


O { X 1U J 




A 7T 


T 1 l TP 

ly + £jp 


O.U ± U.O 


o.yo 


o.y 7 


4.006 




IV p 


T _L TP 
lp + rjy 


1U.Z ± U.o 




y.yU 


1 1 om 

ll.ZUl 




7TV 


i in tp \ 


Z.o ± U.4 


2.72 


2.81 


q ons 
o.zUo 




r 


1 r TP, r\ 

-^[Cy-Ey) 


i a i j- n on 
1.4 ( ± U.zy 


1 AQ 




U. ( ow 




K*°n 
i\ '/ 


1 ( n i TP TP \ 

-y^K^P + Ep — tjy + btiy) 


l.o ± U.4 


1 








K ijj 


1 i i tp \ 
--^(Cy+Ey) 


o o i r\ A 

1.1 ± (J. 4 


Oil 

z.ll 


i on 

1.80 


1.855 




K 4> 


—Ep — SEp 


0.94 ± 0.11 


0.95 


0.90 






K+K*- 


T'y + E' p 


0.20 ± 0.11 


0.25 


0.25 


0.290 




K~K* + 


T'p+E'y 


0.38 ± 0.08 


0.43 


0.43 


0.431 




j\ is. 


E'y — E'p 


< 0.17 


0.08 


0.16 


0.052 




~K°K*° 


a p — h,y 


< 0.09 


n no 
U.Uo 


U.lb 


0.062 




i 




<v U.14 


0.12 


0.12 


n i 






1 l C _L TP IIP i A Q TP \ 
— ^7g V^P + + <£&V + 4<3.E/y j 


/ n in 
< U.1U 


U.UU^i 


n nn^ 

u.uuo 




no 
D 




1 /'/^ , ' OQTP 1 J_ QTP 1 \ 

^=(C P - z,bii p + bh v ) 


zoo 
< z.o 


U.UoO 


U.Uo4 






+ — 


/m I i 7— i / \ 

-(T' y + £ p ) 




0.34 


0.35 


0.485 




7T p + 


-{T'p + E'y) 




0.DZ 


0.61 


0.706 




7T°P° 
r 


1 //"I / 1 (-11 Tp TP \ 
2^ p + L- y — &P — tiy) 




n 1 q 


U. -LU 


n oi a 
u.zlo 






MC, r - C' n + E' n 4- EL 4- 2SE' r ,) 

2\ \f P P V -P / 




0.020 


0.003 


0.013 




rjuj 


— i (C" P + 2C ' v + S'.B ' T/ + 4S , E ' p ) 




0.13 


0.10 






rj'uj 


A-(C7' P - C' v + 4S , J B' V - 2SE' P ) 




0.0007 


0.0003 






rjp° 


— r=f2C t/ — Cp — SE\r) 




0.0039 


0.0015 






rj'p 


i (£" + (7' + 4S'EM 

2 Aj K ^ w P l 




0.012 


0.009 


0.039 




K*+ir- 


"(2? + ^) 




0.029 


0.029 


0.025 




K + p- 


-W + ^p) 




0.016 


0.016 


0.004 




K*°n° 






0.0052 


0.0064 


0.008 




K°p° 






0.0069 


0.0059 






K*°r] 


-^(Cp-^p + ^ + 5^) 




0.0030 


0.0041 






K*°ri' 


^(C P + 2^ + ^ + 45^) 




0.0 


0.0 






K°CO 


^ 2 (Cy + E'p) 




0.0076 


0.0056 


0.002 






E'y + S£p 




0.0 


0.0006 
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TABLE V: (continued). 



Decay Modes 


Representation 


Experimental 


Present ^(xlO -2 ) 


LP [40J 






^(xlO' 2 ) 


FIT A 


FIT B 


i3(xl0" 2 ) 




K 7T+ 


T v + Cp 


1.92 ±0.19 


1.96 


1.96 


1.996 




4- J. 


C p S A p 


0.61 ±0.06 


0.64 


0.62 


0.619 




A p 


Tp + C v 


6.6 ± 2.5 


7.56 


8.43 


12.198 




TT + P° 


1 (rpl , (~t 1 Al A 1 \ 


U.104 ± O.Olo 


u.uoo 


u.uoo 


n i r\A 
0.104 




K+TT*° 
A A 


rp f Al 

1 V A v 


0.42 ± 0.05 


n AA 


n AA 


n A oa 
0.436 




~K° K* + 


rpl Al 
P P 


o. i in i .^i 


1 A Q 
1.4o 


L.ZO 


1.010 




K+p° 

r 


l (n" aii \ 


n 09^ -4- n m 9 


0.030 


0.025 


A A9Q 




K*°ir+ 


(nn i A" \ 

-{Up + Ay) 


n n^R x n nifi 

U.UoD It U.U1D 


094 


0.022 


A 097 

u.uz / 


D+ 


K+<j) 


-(Mr + SA'L) 


< 0.013 


0.0066 


0.0067 






TT + LO 


-4(T' +Cp + A'+A'+ 2SA'p) 




0.57 


0.58 






7]p + 


MT'p + 2C'y+ A'y ± A'p + SA'y) 

-/Q v Jr V V Jr V / 




0.24 


0.43 






v'p + 


— i (T'p - Ct, + At, + A' P + 4SA' V ) 




0.15 


0.15 






TT°p+ 
r 


— J-(T' A-C' a. a' — A' \ 

^{1 p -f- U y "I" /i p /iyj 




0.28 


0.35 


451 

\J. L ±KJ J- 




A p 






n 09^ 


D 099 


042 




7T°K* + 


Lfr"' _ 4'M 




U.Uo 1 


U.UoD 


057 




/ \ UJ 


--^(Cy+Ap) 




n m 9 


n ni 1 

U.Ull 






A Tj 


7§Up - a p + A v + bA v) 




0.015 


0.015 






A r\ 


1 /mil tr\ A II A II A m A II \ 

~M T p + 2A p + A v + ASA v) 




0.00014 


O AAA 1 c 

0.00016 








Cp + Ay 


3.3 ±0.9 


3.34 


3.42 


4.812 




A A 


C V + Ap 


4.3 ± 1.4 


4.98 


4.66 


2.467 




" P 


1 / A A \ 

T2 {Ay-A P ) 


f\ A/lt ( f\ f\l~7\ 

0.06 + (< 0.07) 


U.UO 


U.Uo 






7T (p 


Ty + SAp 


3.6 ± 0.9 


Q HQ 


z.yo 


4.552 




7T A 


(rpl Al\ 


n rc; -X- n 98 


U.oo 


n or 
U.OO 






K+p° 


1 (r" j- a M 


^ n 9Q 


n 1 9 


n 19 

U. 1Z 


u.iyo 


Dt 


K + <j) 


t-i/ \ n i a_ a * i c /t / 

1 y -|- p T /ij/ T O/ip 


< 0.05 


0.032 


0.033 


0.008 




K+uj 


— ^(C"p — A'p — 2SA p ) 




0.40 


0.39 


0.178 




K°p+ 


" -(Tp-Ap) 




0.91 


0.77 


1.288 




7T°K*+ 


-^(CV + A' V ) 




0.13 


0.13 


0.076 




T]K* + 


-y^(Tp + 2C y + A'p — Ay — SAy) 




0.038 


0.047 


0.146 




ri'K*+ 


^(2T'p + C'y + 2A'p + A'y + 4SA'y) 




0.068 


0.059 






K*°K+ 


-(Ty + Cp) 




0.0015 


0.0015 


0.006 




K* + K° 


-(T£ + Cp) 




0.0076 


0.0085 


0.018 



■f The central value of the E791 experiment [39]. 
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TABLE VI: SU(3) flavor symmetry relations of D — > PP decay modes and breaking of the 
relations. A = \V* V us /V* V ud \ « 0.226. k = \V* S V US /V* L V US \ « 4.446. 



SU(3) Symmetry Relations 


LHS of Relations 


FIT a 


FIT f3 


|^(D"^7r+7r-) + v / 2^(I3 ^7r (, 7r )| _ i 
Iv^-Ap+^TT+Tr )! 


1.00 


1.00 


|^(D ^A—7r+)+v / 2^(D ^i? U 7r )l _ ]_ 
|_4(D+-»i?%+)| 


1.00 


1.00 


\\A{D+^TT+Tt') + K,A{D+^K°TT+)\ _ , 

| v / 2K^( J D+^A'+7r0)| 


0.49 


0.79 


\\A(D+^T?'tt+)+V2kA(D+^K+it q )\ 1 
\kA(D+^K u tt+)\ 1 


1 56 


1.11 


|v / 2K^(D+^7r°A'+)-K^(D+^A" ) 7r+)| _ ^ 
\\A(D+^Tt\+)\ 


2.21 


1.82 






T - 1 


1.27 


1.24 




K^(D u ^A-+7r-) 

|^(D U ^7T+7T-)| 


- 1 


1.43 


1.43 


\kA{D"->K+-k-)\ 
|A^(D u ^A--7r+)| L 


1.20 


1.24 


|X4(D°-0?V)| 1 
|^(D"^7r°7r ( ')| ~~ 1 


1.26 


1.12 




\4(£>°-0?V) 


- 1 


1.78 


1.10 


|^(D+^A-+A-")| _ 1 


0.89 


0.86 


|v / 2K.4(£ |+ ^A-+7r )| 
|^(D+-»i^7r+)| 


1.24 


1.24 


|A^l(D+^7f J X+)| _ -, 
\V2A(D+^K+tt°)\ 


1.34 


0.98 


\kA{D+^K u tt+)\ _ 1 


1.08 


1.14 


\\A{D+^Tt\+)\ _ 1 
|v / Z4(D+^7r°7r+)| 


0.55 


0.67 
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TABLE VII: SU(3) flavor symmetry relations of D — > PV decays and breaking of the relations. 
A = \V* s V us /V* s V ud \ « 0.226. k = \V* s V us /V* d V us \ « 4.446. 



SU(3) Symmetry Relations 


LHS of Relations 


FIT A 


FIT B 


\A(D°^tt+K*-)+V2A(D"^tt"K* 1 ')\ _ 1 
\A(D+^tt+K* U )\ 


1.00 


1.00 


\A(d"^p+k-)+V2A(d"^p°k")\ _ 1 

UCD+^"ff"o+1l 


1.00 


1.00 


|.4(D ^i?V)-.4(D+^7r+</>)| -, 
\A(D u ^n+ K*~)\ 


1.00 


1.00 


\A(D°~*7T+K*-)+A(D°^K u <p)\ _ 1 

|.4(D+^7T+</>)| 


0.99 


0.99 


|^(D u ^7r+X*-)-^(D+^7r+</»)| _ 1 
\A(D°^K°<t>)\ 


1.00 


1.00 


\XV2A(Dt~*n+p°)+V2A(D+^n+p°)\ _ j 
|A^(D+^tt+'K* )| 


0.88 


0.88 


|A^(D+^7r+X* U )+v / 2^(D+^7r+p")| -, 
|AV2.4(£>+->7r+p°)| 


0.60 


0.59 


|A^(D+-»7r+i<: )+Av / 2^(D+^7r+ / 9°)| -. 
|v / Z4(D+->ttV , )I 


1.10 


1.10 


|Av / 2^(Dj^ 7 r+p )-v / 2^(D+^7r p+)| _ 1 
|A^(D+-*p+i?°)| 


1.03 


1.03 


\\A(D+->p+K")+V2A(D+-nr°p+)\ , 
\\V2A(D+-+ir+p°)\ 


1 AO 

1.48 


i 1 n 

1.10 


\xV2A(Dt^TT+p°)\-XA(D+^p+K u ) _ 1 
|v / 2^(D+^ 7 r°p+)l 


u.yo 


n q.7 
u.y i 


I i *l) i i ill 

\\A(D+->K+K ")+A(D+->K+K °)| _ 1 
|A.4(D+^+iT U )| 


U.Oo 


n ^7 


|Ayl(D+^7r+_ft:* )-_4(-D + -Kf£T+.F£'* )l _ j 
|XA(D+-^+7T u )| 


117 

i.i i 


i i £ 

1.10 


\\A(D+->K+K )-\A(D+-nr+K )\ _ ^ 
\A(D+^K+K*°)\ 


0.96 


0.97 


\xa(d+^k 1 'k*+)+a(d+^k { 'k*+)\ _ 1 

\XA(D+-+p+K°)\ 


1.09 


1.19 


\XA(D+ ^p+Tt ')-A(D+ ^T?' 'K*+)\ _ 1 
\\A(D+^K u K*+)\ 


1.01 


0.97 


\XA{D+ ^T?' K*+)-XA(D+ ^p+ K")\ _ 1 
\A(D+^K i, K*+)\ 


0.98 


0.96 


\XA(D+^K i 'K*+)+y/2A(Df^K*+TT°)\ _ -. 

IA\/2.4.CD + — »7r+o 1l ~ 


1.42 


1.31 


IAv^DJ^ttVO+v^DJ-^'+tt )! r 

, rt — 1 

|A^(D+-»K 


0.92 


0.95 


|Av / 2^(D+^7r+p°)+A^(D+^if , ^*+)| -, 
\V2A(D+^K*+n°)\ 


1.11 


1.07 


\A(D+-> K°p+)\ , 
|.4(D+^i?V*+)| 


n qi 


u.oy 


\A(D+^n+K*U)\ _ 1 
|^(D+^AT+K* U )| 


0.93 


0.94 


|^(D ,J ^K+^*-)| -, 
M(D"-»A-+p-)| i 


1.05 


1.05 


|A^(D u ^7r+K*-)| n 
\A(D"^ir+p-)\ L 


1.05 


1.05 


\A{D"^K + K—)\ -, 
|A^(D u ^7r+E"*-)| — L 


1.14 


1.14 


\A{D u ^K-K*+)\ n 
\kA(D u ^tt-K*+)\ - 1 


1.08 


1.08 


\XA{D»^K-p+)\ _ 1 
|_4(D'J->7r-p+)| _ X 


1.09 


1.09 


|^(D U ^A"-E"*+)| -, 
\XA(D"^K-p+)\ - 1 


1.08 


1.08 



